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Abstract 


A  recent  result  by  Panjer  provides  a  recursive  algorithm  for  the  com¬ 
pound  distribution  of  aggregate  claims  when  the  counting  law  belongs 
to  a  special  recursive  family.  In  the  present  paper  we  first  give  a 
characterization  of  this  recursive  family,  then  describe  some  general¬ 
izations  of  Panjer 's  result. 


1.  Introduction 


Let  u  be  the  Lehesgue  or  the  counting  measure  on  (0,®)  ,  and 
let  x1#x2, . . .  be  independent,  identically  distributed  random 
variables  (the  independent  severities)  with  cumulative  distri¬ 
bution  F  and  generalized  density  f  : 

F (x)  =  /  f(y)  du(y)  . 

(0  ,x  ] 

Let  ft  be  a  random  variable  (the  claim  number) ,  independent  of 
the  x^s  ,  defined  on  the  non-negative  integers  with  probabili¬ 
ties  : 


pn  =  Pr (n  -  n)  . 

Then  the  generalized  density  g  of  the  random  sum  (the  aggregate 
claims) 


n 

3*1 

i=l 


(wetacitly  assume  s  is  zero  if  n  is) 
has  an  atom 

(1.1)  g (0)  *  pQ 

at  zero,  and  for  s  >  0  the  form 

ao 

(1.2)  g(s)  *  :  pn  fn*(s)  , 

n=l 

where  f  denotes  the  n-th  convolution  of  f  .  This  formula 
is  extremely  difficult  to  compute  because  of  the  high-order  con¬ 
volutions;  only  a  few  closed-form  solutions  are  known. 

1 


Panjer  (1981)  has  shown  that,  if  there  exist  constants  a  and 
b  such  that 


(1.3) 


pn  =  pn-l  (a  +  n}/ 


(n  —  1,2,...) 


then 


(1.4) 


g  (s) 


p.  f  (s)  +  /  (a  +  b  f)  f  (x) 

(0 ,  s ) 


g(s-x)  du (x)  . 
(s  >  0) 


The  importance  of  this  result  is  that,  when  f  is  discrete,  the 
successive  values  of  g  can  be  recursively  calculated.  We  now 
consider  various  aspects  of  the  relation  between  the  recursions 
(1.3)  and  (1.4),  and  then  provide  a  variety  of  generalizations. 


2.  Characterization  of  the  counting  distribution 


The  following  theorem  characterizes  the  class  of  counting  densi¬ 
ties  pn  satisfying  (1.3);  it  is  essentially  given  in  Johnson  & 
Kotz  (1969). 

Theorem  1 

Assume  that  (1.3)  holds.  Then  we  must  have  one  of  the  four  cases: 


(  0  (n  =  0 ) 
(n  >  0) 


(2.2) 


(X  >  0) 
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(2.3) 


P„  =  p"<l-p>“ 


(a  >  0,  0  <  p  <  1 ) 


(2.4) 


pn  =  (n5 


n , ,  . N-n 

P  (1-P) 


(0  <  p  <  1,  N  =  1,2,...) 


Proof .  To  avoid  negative  probabilities  we  must  have  a  +  b  i  0  . 

For  a  +  b  =  0  we  get  the  degenerate  case  (2.1).  For  the  rest 
of  the  proof  we  assume  a  +  b  >  0  .  If  a=0  we  get  the 
Poisson  density  (2.2)  with  X  =  b  .  For  a  >  0  we  introduce 
a  =  (a+b)/a  and  get  from  (1.3) 


,3-t-n-L  n 
pn  '  po  1  „  >  a 


In  order  that  E  ,  p  <  1  ,  we  must  have  a  <  1  .  Then  we  get 
n=l  n 

the  negative  binomial  (Pascal)  density  (2.3)  with  p  =  a  . 

Finally,  assume  a  <  0  .  Then,  to  avoid  negative  probabilities, 
there  must  exist  a  positive  integer  N  such  that 
a  +  b/(N+l)  *  0  ,  that  is,  N  =  -(a+b)/a  .  With  p  =  -a/(l-a) 
we  get  the  binomial  density  (2.4). 


We  have  now  proved  the  theorem. 


Q .  E .  D  . 


The  allowed  regions  for  (a,b)  are  illustrated  in  Figure  1  , 
which  is  inspired  by  Johnson  &  Kotz  (1969,  p.  42). 

Remark .  For  the  case  a  <  0  Johnson  4  Kotz  (1969,  p.  41)  also 
develop  a  distribution  for  the  case  when  -(a+b)/a  is  not  an 
integer,  by  letting  Pn  =  0  when  a  +  (b/n)  <  0  .  However,  that 
distribution  does  not  satisfy  (1.4)  as  we  then  must  have  that 
(1.3)  holds  for  all  n  >  0  .  A  modified  version  of  (1.4)  allow¬ 
ing  such  "generalized  binomial"  distributions  will  be  given  in 
Section  5.  However,  this  version  seems  in  most  cases  to  be  more 
complicated  than  direct  computation  of  (1.2).  For  the  binomial 


BINOMIAL 


b 


Figure  1.  Permitted  (a,b)  parameter  space 

for  recursion  (1.3)  (The  dotted  area 
denotes  the  increase  obtained  by 
recursion  (3 .10)  .) 
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\j 

distribution  we  have  that  Pr  (n  >  N)  =0  ,  but  as  (*  )  =  0  for 

n 

n  >  M  we  can  let  be  defined  by  (2.4)  for  all  the  non-nega¬ 

tive  integers. 


3.  Generalizations 


We  first  introduce  some  notation:  if  are  given 

quantities,  then  we  let 


z 


nZ 


z  . 
i 


denote  the  sum  of  the  first  n  elements. 

Assume  there  exists  a  function  h  :  {(x,s)  :  0  <  x  <  s }  -  H  , 

satisfying  the  condition  that 


(3.1) 


~i  (h  (x.  ,s)  ]  x  _  *  s)  =  in  (n  -  2,3,...) 

l  nl  n 

are  independent  of  s  . 


Then  we  have  the  following  generalization  of  Panjer's  result: 


Theorem  2 
If 

(3.2)  Pn  =  Pn_j_  mn  '  (n  =  2,3,...) 

with  the  sequence  satisfying  (3.1),  then 

(3.3)  g(s)  =  p  f(s)  +  /h(x,s)  f(x)  g(s-x)  du (x)  .  (s  >  0) 

(0,s) 


Proof:  We  have  for  s  >  0 


-o- 


g  (s) 


(<5) 


n*L 


p.  f(s)*r”,p  .  m  fn(s)  = 

r 1  n*2  n-1  n 

p  f(s)  P„  ,  h(x,s)  f  ( x )  f<n  (s-x)  du(x)  = 

1  n"2  n_1 ( 0 , s ) 

f(x)ri“  .  p  fn  (S-x)]  du (x)  = 
n  =  l  n 

f(x)  g(s-x)  du(x)  . 

Q.E.D. 


p .  f  ( s )  ♦  /  h  ( x ,  s  ) 

( 0  ,  s ) 

p  f  (s )  +  /  h (x,s  ) 
(0 ,  s ) 


It  is  clear  that  if  the  functions  h^  and 
(3.1),  then  for  all  constants  c^  and  c^ 
c^  h^  +  c2  h2  satisfies  (3.1). 


h2  both  satisfy 
the  function 


For  all  constants  a  and  b  we  clearly  have 

x,  . 

(3.4)  S(a  +b  x  -s)  =  a  +  2  ,  (  n  =  2,3 - ) 

s  t\L  n 


independent  of  s  .  Hence  the  kernel  in  (1.3)  , 


h (x,s)  =  a  +  b  -  , 


is  a  special  case  of  (3.1)  with 


(3.5) 


.  b 

m  =  a  +  - 
n  n 


(n  =  2,3,...) 


The  following  example  gives  a  distribution  satisfying  (3.1)  with 


m  satisfying  (3.5),  but  not  covered  by  Panjer’s  result, 
n 


Example  1.  Consider  the  logarithmic  counting  density 


(3.6) 


pn  = 


(n  =  0) 


n 

E_ 


(0  <  p  <  1) 


l«.n(l-p)  |  n 


(n  =  1,2,  .  .  .)  . 


Then  we  have 


? 


n 


?n-l 


(n  =  2,3,  .  .  .)  , 


that  is,  m  =  p [  1  -  (1/n)  ]  ;  a  =  p  ;  b  =*  -p  ,  and  for  s  >  0 
n 

we  obtain 


g(s)  -  p.  f (s )  +  P  /  (1  -  f )  f(x)  g(s-x)  du(x)  . 

1  (0,s) 

The  difference  from  Panjer's  result  is  that  (1.3)  does  not  hold 
for  n  =  1  . 


Theorem  3 

Assume  that  (3.1)  is  satisfied  for  the  distributions  given  by 

Pr(x.=l)  =  1  -  Pr(xi=2)  =  p  ,  (pc  [ 0 , 1 ! )  . 

Then  there  must  exist  constants  a  and  b  such  that  (3.5)  is 
satisfied . 


Proof .  For  p  =  1  and  p  =  0  we  get 


(3. 

.7) 

m  =  h (1 ,n ) 
n 

(3 

.3) 

m  =  h (2 , 2n) 
n 

respectively. 

Assume  p  e  (0,1)  ;  u  =  2,3,...  ;  n  =  u  ,  u  +  1  ,  ...  ,  2u  . 

Then 


-s- 


fn-1  \ 

\ 2u-y-n+l ) 


(3  9)  f(v)  f(n  L>  (2u-y)  _  V  2u-y- 

’  ( 2u-n) 


(y  =  1,2) 


By  using  (3.7),  (3.8),  and  (3.9)  in 


f  (1 )  f(n-1,*(2u-l)  ...  „ 

m  =  - — -  h(l,2u) 

n  -n* 


f  (2u) 

(n-1 ) * 

f ( 2 )  ftn  (2u-2 ) 

n  * 

f  ( 2u ) 


h (2 , 2u) 


we  obtain 


m  =  a  +  —  i 
nun 


with 


a  =  2m-  -  m  , 

u  2u  u 


bu  =  2u(mu':n2u)  • 


As 


.  u+l 

mu+l  au+l  u+l  au 


u+l 


mu+2  “  au+l  '  u+2 


u+l 

+  — — —  =  a 


b 

_ u_ 

u+2  ' 


we  must  have  a  ..  =  a  and  b  =  b  for  all  u  ,  Lnat  is, 
u+l  u  u+l  u 

there  exist  constants  a  and  b  such  that  (3.5)  is  satisfied. 

Q .  E .  0 . 


Theorem  3  says  that  if  (3.1)  is  to  hold  for  a  class  of  two-point 
distributions  F  the  sequence  (mn }  must  satisfy  (3.5).  This 
result  clearly  implies  tnat  if  '3.1)  is  to  hold  for  all  distri¬ 
butions  on  (O,00),  the  sequence  {mn}  must  satisfy  (3.5).  Be¬ 
cause  of  this  fact  we  restate  Theorem  2  for  this  particular  class 
of  counting  distributions. 
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Theorem  2 ' 


If 

(3.10)  Pn  =  pn-l^a+n^/  *n  *  2,3,...) 

then  for  all  severity  distributions  F  we  have 


(3.11)  g(s)  =  p,  f  ( s )  +  /  (a+b£)  f(x)  g(s-x)  du(x)  .  (s  >  0) 

(0 ,  s ) 

We  close  this  section  by  comparing  the  class  of  counting  distri¬ 
butions  defined  by  (1.3)  (that  is,  the  class  given  in  Theorem  1) 
to  the  class  defined  by  (3.10).  Clearly  the  latter  class  con¬ 
tains  the  former  one.  As  in  the  latter  class  the  recursion  may 
start  at  one,  the  restriction  a  +  b  i  0  may  for  a  >  0  be  re¬ 
placed  by  the  weaker  condition  a  +  b/2  *  0  .  Hence,  the  permit¬ 
ted  parameter  space  is  now  increased  by  the  dotted  region  of 
Figure  1. 

As  pQ  may  now  be  chosen  (relatively)  freely,  the  counting 
distribution  is  no  longer  uniquely  determined  by  (a,b)  .  For 

(a,b)  being  in  the  permitted  region  for  recursion  (1.3),  ex¬ 
cluding  the  line  a  +  b  =  0  ,  the  permitted  class  consists  of 
the  distributions  given  by 


(3.12) 


r  p  +  (l-p)n 

(n  =  0) 

pn  =  * 

l  . 

(n  =  1,2,  . 

where  {  tt  } 
n 

0  is  chosen 
fies  (3.10) 
case  (3.11) 


is  a  counting  distribution  satisfying  (1.3),  and 

such  that  0  «  1  and  p_  i  0  .  p  deary  satis- 

o  n 

with  the  same  (a,b)  as  for  n  .  in  the  discrete 

n 

may  under  the  present  conditions  be  written  as 


*  (a+b)  p  f (s) 


+  Z  (a  +b  £  ) 
x*  1  s 


f(x)  g(s-x)  . 


g  (s) 


(s  >  0) 


For  a  +  b  = 
is  given  by 


Pn  given  by  (3.6).  A  counting  distribution 


0  the  permitted  class  of  counting  distributions 

(3.12),  with  the  obvious  restrictions  on  p  ,  and 

p  }  of  the  form 
n 

(3.12)  may  be  interpreted  as  a  weighted  (in  a  general  sense, 

as  p  may  be  negative)  distribution  of  the  distribution  { tt  } 

n 

and  a  distribution  concentrated  at  zero.  Then  the  aggregate 
claims  distribution  must  be  the  analogous  weighted  distribution 
of  aggregate  claims  distributions,  and  if  the  aggregate  claims 


distribution 


corresponding  to  rr  is  known,  we  may  find 


the  aggregate  claims  distribution  g  corresponding  to  p  by 


gp(s)  = 


ro  +  ( 
l  (l-o) 


(l-o) 


v°> 


gT(s) 


(s  =  0) 
(s  >  0) 


4.  Results  on  specific  severity  distributions 


From  (3.4)  and  Theorem  3  we  see  that  if  (3.1)  is  going  to  be 
satisfied  for  all  F  ,  then  the  sequence  {nO  must  satisfy 
(3.5).  However,  for  specific  classes  of  F  there  exist  other  m 
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The  following  obvious  result  is  interesting  in  this  connection. 

Theorem  4 

Let  v  be  a  function  such  that  v(j^,xnr)  is  independent  of 

x  _  for  all  n  .  Then  (3.1!  holds  for  any  h  that  can  be  written 
nE _ _  _ 

h(x,s)  =  k  (v  (x,s))  with  £(h  (xx  ,xn  )  )  existing  for  n  =  2,3,... 

Example  2.  Assume  that  x1,x2,...  are  gamma -distributed  with 
parameters  (a,v)  .  Then  x^/x^  is  independent  of  x^  and 
beta -distributed  with  parameters  (v,(n-l)v)  .  Hence,  by  Theorem 

4,  all  h  (x,  s )  =  k(x/s)  with  (x,/xnr )  )  existing  for  all  n 

satisfy  (3.1).  In  particular,  if 

k(z)  =  zu(l-2)v 


we  get 


T  (nv)  r  (v*u?  F  (  jn-l)  v+v) 
mn  F((n-1)  v)  F(v)  F (nv+u+v) 

For  v  =  0  and  u  positive  integer  this  gives 

u-1  u-1  a. 

0.  -  -1  =  I  -±r 

n  i-o  nv  1  i-0  nv+1 

for  some  aQ,...,au_^  independent  of  n  .  Hence,  for  any  positive 
integer  u  there  exist  constants  c^,...,c  ^  such  that 

U+1  i 

k (z)  =  Z  c.  zx 

i-1  1 


gives 


m 


n 


1 

nv+u 
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Example  3 .  Assume  that  the  counting  density  is  hypergeometric 


m  I 
n/\ N-n  / 


where  the  positive  integer  parameters  (m,M,N)  satisfy  N  <  M  ; 
m  s  M  -  N  .  For  n  >  0  we  have 

_  (m-n+I ) (N-n+1) 

^n  ^n-l  n(M-m-N+n)  ' 

which  may  be  written 

^n  ^n-1  +  n  n+M-m-N^  ' 


with 


a  =  1  , 

.  _  (m+l) (N+l) 

N-M+m 

(M-m+1) (N-M-l) 
c  N-M+m 

Now,  assume  that  the  x^s  are  gamma -distributed  with  parameters 
(a,v)  ,  where  v  is  a  positive  integer.  As  we  may  write 

c  _  _ cv _ 

n+M-m-N  ~  nv+(M-m-N)v  ' 

by  Example  2  we  can  find  a  function  h  such  that  Theorem  2  is 
satisfied. 

The  extension  to  the  eccentric  hypergeometric  distribution  (see 
Sverdrup  (1976),  with  counting  density 


I 
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(X  >  0) 


where  pn  is  given  by  (4.1),  is  obvious. 

Similar  approaches  are  possible  for  the  following  counting 
distributions,  described  in  Johnson  &  Kotz  (1969):  the  displaced 
Poisson  distribution  (p.  113);  and  the  Yule  distribution  with 
generalizations  (pp.  244-251). 


5.  Recursion  on  a  limited  range 


In  the  previous  cases  we  have  assumed  that  the  o  can  be  com- 

**  n 

puted  recursively  for  n  >  1  .  The  following  Theorem  5  extends 
this  to  the  case  when  the  recursion  holds  only  for  n  >  K  with 
K  s  1  . 

Let 


gK(s)  = 


Pn  «  (S) 


n=K 


Then 

K-l  „ 

g(s)  =  I  Pn  fn  (s)  +  g  <s)  . 
n=0 


Theorem  5 


Assume  that 


P_  *  P„  ,  m  .  (n  *  K+l , K+2 , . . . ) 

n  n- 1  n 

with  mn  given  as  in  (3.1)  .  Then 
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(5.1) 


gK(s) 


p  fK  (s)  +  f  h(x,s)  f(x)  g  (s-x)  d;i(x)  . 
K  (0 ,  s )  K 


(The  proof  goes  as  in  Theorem  2  and  is  omitted.) 

The  difference  from  the  underlying  assumptions  of  Theorem  2  is 
that  (3.1)  and  (3.2)  do  not  need  to  hold  for  n  r  X  .  If  (3.1) 
holds  for  all  n  i  2  ,  insertion  of 


gR(s)  =  g(s) 


K-l 

£ 

n=*l 


-n*  ,  . 
f  (s) 


(s  >  0) 


in  (5.1)  gives  the  final  recursion: 


(5.2) 


g  (s)  =  p,  f  (s) 


n=2 


(P  “D 
r  n  “  n  — 


n-1 


V 


(s) 


+  f  h(x,s)  f  (x)  g (s-x)  du(x)  .  (s  >  0) 

(0 ,  s ) 

(The  summation  is  zero  if  K  *  1  . )  Compared  to  (3.3)  we  have 
now  got  the  summation  as  a  correction  term,  since  this  would  be 
zero  if  Pn  -  ?n-i  mn  =  0  for  n  =  2,...,K  . 


For  the  special  case  of  Theorem  5  with  pQ  =  =  ...  =  px_^  =  0 

(truncation  from  below)  g  (s)  *  g(s)  ,  and  (5.1)  gives 


(5.3)  g (s)  *  p„  fK  (s)  +  f  h(x,s)  f(x)  g(s-x)  du(x)  . 

(0,s) 
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We  shall  now  see  what  happens  if  the  counting  distribution  is 
truncated  from  above.  Assume 

f  0  (n  =  0 , . . . , K-l ) 

Pn  *  1  Pn-1  mn  (n  =  K+1,...,L) 

i 

^  0  (n  =  L+l , . . . ) 

Then  for  s  >  0  we  get 

(5.4)  g (s )  =  pK  fK  (s)  -  PL  mL+1  f(L+1)  (S) 

+  /  h(x,s)  f(x)  g(s-x)  du  (x)  . 

(0 ,  s ) 

Unfortunately, in  this  formula  we  need  high-order  convolutions  of 
f  .  These  can  be  rather  complicated  to  compute,  except  for  some 
cases  where  we  have  simple  closed-form  expressions  (gamma, 
Poisson,  binomial,  negative  binomial  distributions).  In  some 

cases  the  factor  p  mT  .  makes  the  correction  term  negligible. 

^  ^  ( L+l )  * 

Another  possibility  is  for  large  L  to  approximate  f  by 

a  (possibly  discretized)  normal  density.  Otherwise  it  is  probably 

more  efficient  to  compute  g  from  the  basic  definition  (1.1). 


6.  Extension  to  non-positive  discrete  values 


We  now  leave  the  assumption  that  the  x^  s  are  distributed  on 
((),<*■)  and  assume  that  they  are  distributed  on  the  set  of  all 
integers : 


f (x)  =  Pr(x=x)  .  (x  =  . .  .-2,-l,0, 1 ,2, . . . ) 


Then  (1.1)  must  be  replaced  by 


We  further  assume  that  the  counting  distribution  satisfies  the 
recursion  (1.3),  and  analogously  to  Theorem  2  we  obtain 


(6.2)  s  g(s)  =  Z  (as+bx)  f(x)  g(s-x)  . 

X--00 

If  x,  only  takes  on  zero  plus  positive  values,  so  does  s  ; 
then  f  (0)  =  [ f ( 0 ) ]  ,  and  the  sum  in  (6.1)  can  be  carried  out 

explicitly  (see  the  probability  generating  function  for  the  count¬ 
ing  distribution  in  Johnson  4  Kotz  (1969)).  We  then  get  the  re¬ 
cursive  system 


(6.3) 


U 


9(0) 


r 

l-a  f(0) 
1  -a 


i 


.  -b[l-f  (0)  ] . 

e  , 

V 


(a  *  0) 

(a  *  0) 


g  (s) 


l-a  f(0) 


Z  (a  +  b  f (x)  g (s-x)  . 

x=l 


(s  >  1) 


The  case  where  the  x^  can  take  on  negative  values  is  difficult 
because  one  cannot,  in  general,  find  suitable  starting  values 
for  s  in  (6.2). 

However, in  the  case  where  pn  is  Poisson  with  parameter  X  (2.2) 
the  density  g  can  be  computed  by  two  applications  of  (6.3)  plus 
a  convolution.  Let 

k*  =*  max(0,x^) 

(6.4)  (1-1,2,...) 

SL  =  max  (0  ,  -x^ )  , 
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and  we  have 

-+  5  .+  5  . 

s  =  E  x .  ;  s  =  E  x .  < 

i=l  1  i=l  1 

(6.5) 

s  *  s  -  s 

to  US 

AndrS  Dubey  has  pointed  outVthat  when  n  is  Poisson  distributed, 

then  §+  and  s  are  independent.  Let  x+  and  x.  have 

+  —  1  +  1  _ 

densities  f  and  f  ,  respectively,  and  s  and  s  have 

densities  g+  and  g  ,  respectively.  Then  g+  and  g~  are 

computed  independently,  using  (6.3),  with  a  =  0  ,  b  =  >  ,  and 

the  corresponding  f+  or  f  Then  g  for  the  total  sum  is 

computed  by  the  convolution 

oo 

(6.6)  g(s)  =  E  g+(x)  g” (x-s ) 

x=max(0 ,s) 

(6.2)  can,  in  principle,  also  be  solved  for  p^  binomial,  if 
f(x)  is  defined  over  (-K, -K+l , . . . )  ,  for  in  that  case  there  is 
a  largest  negative  value  of  the  sum,  s  =  -NK  ,  and  (6.2)  can 
be  rearranged  into  a  true  recursive  form. 

Remembering  that  p  =  -a/(l-a)  and  N  *  -(a+b)/a  ,  we  get 
the  recursive  system: 


(s  >  -NK) 
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Of  course,  if  K  is  very  large,  there  are  obvious  problems  with 
round-off  error  accumulation,  especially  if  f(-K)  and  the  nearby 
values  are  very  small.  We  remind  the  reader  that  this  problem  can 
occur  with  any  recursive  scheme  described  in  this  paper  where  the 
range  of  discrete  severities  is  large. 

There  remains  the  case  of  pn  negative  binomial  (2.3)  for  which 
it  does  not  seem  possible  to  give  a  simple  procedure  for  negative 
55!  s  .  Of  course,  in  this  and  in  the  other  cases,  one  can  think 
of  various  iterative  schemes  for  (6.2)  which  would  converge  to 
the  correct  density. 
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Abstract 


A  traditional  actuarial  method  for  the  difficult  task  of 
finding  the  exact  distribution  of  a  heterogeneous  portfolio 
approximates  the  distribution  with  a  compound  Poisson  law 
with  identically  distributed  risks.  This  paper  shows  that  a 
Binomial  compound  law  provides  a  better  match  to  the  second 
moment  of  the  distribution,  thus  giving  a  better  approximation 
while  retaining  a  simple,  recursive  algorithm  for  calculating 
the  distribution.  A  modified  Binomial  compound  law  further 
refines  the  approximation,  with  slight  additional  work. 


Introduction 


Let  (x^;i=l,2, . . .N)  be  a  fixed  set  of  independent,  non-identically 
distributed,  integer-valued  random  variables  for  which  the  probability 
that  any  x^O  is  significant;  we  wish  to  find  the  distribution  of  the 
sum  y=x^+X2+. • *+x^.  In  principle,  the  discrete  density  of  y  is  calcu¬ 
lated  as  the  N-fold  convolution  of  the  discrete  densities  of  the 
individual  x^;  however,  this  task  is  already  very  time-consuming  on 
digital  computers  for  N  larger  than,  say,  1,000,  if  the  x^  take  on 
more  than  a  few  different  values. 

An  approximate  method,  used  for  many  years  by  actuaries,  utilizes 
the  fact  that  many  terms  in  the  sum  may  have  value  zero  and  computes  y 
as  if  it  were  the  sum  of  a  random  number  of  identically  distributed 
random  variables;  in  this  method,  the  first  moment  of  y  is  matched 
exactly,  and  the  second  moment  is  matched  approximately. 

In  this  paper,  we  present  an  improved  approximation  method  that 
provides  a  much  closer  fit  to  the  second  moment,  yet  maintains  a  simple, 
recursive  algorithm  for  computing  the  density  of  the  random  sum. 

Limited  computational  experience  indicates  that  the  approximation  to  the 
distribution  of  y,  and  other  functions  of  y,  are  much  closer  to  their 
true  values  than  in  the  classical  method. 


The  Heterogeneous  Portfolio 


:or  tne  moment,  we  assume  that  the  x^  take  only  r.on-negative 
values  m  the  range  '0,1,...R.‘,  and  we  separate  the  given 
discrete  density  of  as  follows: 


(1) 

Pi 

Pr 1X^=0;  =  1  -  q^ 

t 

( i=l , 2 , . 

.  .N) 

(2) 

fi(x)  = 

?r(xi=x  x^>0;- 

• 

(x=l , 2  ,  . 

.  .R) 

(This  is  a  traditional  notation) . 

We  wisn  to  calculate  the  discrete  density  g  of  the  sum 
N 

(3)  y  =  -  x.  ,  (y=0 ,1, . . .NR) 

1=1  1 

which  is  given  exactly  by  the  N-fold  discrete  convolution: 

N 

(4)  g (y)  =  ^  .'p.  <y)  +  g,  f.  (y)  1  » 

where  :(y)=l  if  y=0,  and  is  zero  otherwise.  The  approximation 
to  be  described  requires  that  "most"  of  the  p^  be  "rather  large". 

Denote  the  first  two  moments  of  the  positive  part  of  the 
random  variables  by 

_  -  -  "> 

(5)  m.  =  Eix.,x.>0)  =  ixf.(x)  ;  v.  =  Vlx.  x.>0;  =  I (x-m, ) “f . (x) 

ill  i  ill  j-i 

Then  it  is  easy  to  show  that  the  first  two  moments  of  the  sum 
y  are: 

N 

(6)  E(y)  =  I  q.m. 

i  =  l  1  1 

N  N  2 

U(y)  =  _  q . v -  +  _  p  q  .  (m.  ) 

1=1  11  1=1  1  1 


(7) 
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The  evaluation  of  g(y)  is  often  required  for  insurance  risk 
portfolios ,  where  i=l,2,...N  indexes  the  policies  m  the  portfolio, 
assumed  independent;  p  is  the  usually  significant  no-claim 
probability  during  a  certain  exposure  period ;  is  the  probability 
of  at  least  one  claim;  and  f^(x)  is  the  density  of  aggregate 
claims  during  the  exposure  period  for  policy  i,  given  that  at 
least  one  claim  occurs. 

The  situation  is  particularly  simple  in  life  insurance,  as 
usually  just  one  claim  occurs  at  death,  and  the  f^(x)  are  often 
only  one-  or  two-point  densities  (e.g.,  the  face  value  of  a 
policy  i  payable  at  the  death  of  the  assured,  who  has  mortality 
rate  in  this  exposure  period) .  Often,  only  the  change  from 
one  exposure  period  to  the  next.  Approximation  methods  have  become 
less  important  in  such  simple  cases,  especially  with  N  small,  as 
modern  computers  can  often  calculate  the  exact  convolution  (4) 
directly.  However,  for  large  portfolios  with  arbitrary  f . (x) ,  the 
problem  of  approximating  g(y)  still  remains.  Most  approaches  have 
been  based  upon  moment-matching,  using  (6)  and  (7) . 


The  Collective  Risk  Model  as  an  Approximation 


One  useful  idea,  from  both  the  theoretical  and  computational 
points  of  view,  is  to  approximate  the  inhomogeneous,  fixed  portfolio 
by  a  homogeneous  risk  collective ,  in  which  we  replace  the  individual 
policies  by  a  mass  of  similar,  anonymous  policies,  and  assume  that 
y  is  the  sum  of  a  random  number ,  n,  of  independent  claims  that  are 
identically  distributed  samples ,  ( w^ ,w2 , . . . w~ )  of  a  positive  random 
variable,  w,  with  some  prototypical  claim  density,  f(w).  If 


(8)  *  =  Prih=n;  (n=0,l,...)  ;  f(w)  *  Pr{w=w;  (w*l,2,...) 

then  this  leads  to  the  well-known  compound  law  of  risk  theory: 


t9) 


(y> 


n=l 


-  '  e  ,  ,  -  n* 

■n- f (y)  i 


g(y) 
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The  rationale  for  this  approximation  is  easily  seen.  If  the 
p,  are  significant,  then  the  sum  y^k^x,  * .  .  .  -t-x^  will  have  a 
varying  number  of  non-zero  terms;  the  sum  could  then  be  repre¬ 
sented  by  y=w1+w2+. . . where  these  all-positive  terms  could 
be  considered  to  be  identically-distributed  samples  from  some 
"representative’’  claim  density,  calculated  by  weighting  each 
fi(x)  according  to  its  probability  of  occurence,  q^. 

If  the  prototypical  claim  moments  are: 


(10) 

m  =  £  ( w ) 

=  Zwf(w)  ;  v  =  i/(w)  =  I(w-m)^f(w) 

/ 

then 

the  moments 

of  the  random  sum  (9)  will  be: 

(11) 

£(y) 

=  £ (n) m  , 

(12) 

i/(y) 

=  £(h)v+l/(n)m2 

For  a  good  approximation,  the  moments  (11)  (12)  must  be  matched 

as  closely  as  possible  to  the  exact  values  (6)  (7),  so  that  g(y) 

and  related  functions  calculated  via  (9)  will  match  values 
calculated  via  (4) . 

We  are,  of  course,  free  in  devising  an  approximation  to 
choose  "  and  f(w)  in  any  way  we  choose.  But  the  most  natural 
way  to  fix  the  prototypical  claim  density,  consistent  with  the 
risk  theory  interpretation,  is  as  the  weighted  sum: 

Iq  f  (w) 

(13)  f(w)  »  — — -  ( w=  1 , 2 _ R) 

-*L 

(This  choice  is  also  invariant  under  pre-aggregation  of  the 
policies  in  a  consistent  way,  for  example,  by  lumping  together 
all  policies  with  the  same  single  face  value  and  adding  their 
q^’s).  With  this  choice,  the  moments  of  w  become: 

(14)  m  =  ( :qimi) / ( Zq^ )  ; 

2  =  +  :qi(mi)  2!/ Cqj) 


(15) 


v  ♦  m 
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If  (11)  and  (12)  are  to  be  matched  exactly  to  (6)  (7),  then 
this  implies  that  the  counting  density,  ^  ,  must  be  chosen  so 
that: 

(16)  E(n)  =  , 

and 

(17)  I /(h)  *  :qj  -  ZqJtnij/m)2 

The  mean  of  n  is  just  the  mean  number  of  positive  terms  in  (3): 
however,  the  variance  of  counts  is  not  the  variance  of  the  number 
of  positive  terms,  Iq^Pj,  unless  the  policies  have  identical 
face  values.  This  is  because  we  are  matching  moments  between  two 
different  models,  one  where  the  sampling  is  without  replacement, 
and  another  where  the  sampling  is  independent. 

Note  that,  in  certain  unusual  cases  where  the  are  small  and  the 
are  quite  different  from  another,  V(n)  in  (17)  may  be  negative;  in  other 
words,  the  approximation  cannot  be  used.  For  instance,  if  N*2,  m^“l, 
m^”?,  and  q^-q^q,  then  we  find  that  q  must  be  smaller  than  0.64  to  obtain 
a  positive  variance.  This  makes  precise  our  earlier  remark  that  most  of 
Che  p^  should  be  rather  large. 


The  Poisson  Counting  Distribution 


A  good  theoretical  case  can  be  made  for  the  Poisson  density: 


(18) 


1 

n 


n! 


(n*0 , 1 , . . . ) 


as  an  appropriate  choice  for  the  counting  law;  Gerber  (1979) 
presents  an  argument  based  on  a  limiting  result  from  the  fixed 
portfolio  model,  as  well  as  an  argument  based  upon  a  dynamic 
portfolio,  in  which  claiming  policies  are  immediately  replaced 
by  equivalent,  non-claiming  policies.  (16)  then  leads  to  the 
natural  choice 


(19) 
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But  from  (17)  it  can  immediately  be  seen  that  the  Poisson 
assumption,  which  means  l/(n)*\,  leads  to  too  large  a  value  of 
i /(h)  for  the  second  moments  (7)  and  (12)  to  match.  In  fact, 
the  collective  approximation  will  now  have  a  variance 

(20)  >/(y)  =  :<3ivi  +  ^i(mi)2 

which  is  greater  than  the  correct  value  (7)  by  the  amount 

-22 
-*i"i  * 

Another,  less  critical,  problem  is  that  the  probability  of 
no  claim  in  the  risk  approximation: 

(21)  g(0)  a  =  e  ^  =  e  , 

is  termwise  greater  than  the  true  value  from  (4) : 

N 

(22)  g (0 )  -  2  p 

j-1  3 


Discussion 


In  addition  to  having  a  good  fit  between  the  approximation 
and  the  original  model,  we  would  also  like  to  have  the  computation 
of  g(y)  via  (9),  and  of  related  functions,  to  be  efficient; 

Gerber  (1980)  describes  some  of  the  traditional  approximations  to 
the  compound  Poisson  law  which  have  been  used  by  actuaries. 

However,  a  simple  recursive  scheme  for  the  Poisson  case, 
apparently  due  to  Adelson  (1966),  has  recently  been  promoted 
as  the  most  efficient  solution  to  (9).  In  our  notation,  it  can 
be  shown  that: 


g(0)  =  e 


(23) 


min(y,R) 

g(y)  *  ( \/y) l  xf(x)g(y-x)  .  (y=l,2,...) 

x=l 

T.iis  enables  exact  values  of  g(0),  g(l)  ,  g(2),  ...  to  be 
calculated  successively,  in  a  number  of  steps  much  less  than 
direct  ways  of  calculating  (9).  A  simple  proof  of  (23),  due 
to  Buhlmann  and  Gerber,  can  be  found  in  Gerber  (1980) . 
Applications  can  be  found  in  Panjer  (1980)  and  in  Held  (1980). 


More  recently,  Panjer  (1981),  has  extended  the  recursive 
computation  of  g(y)  to  a  larger  class  of  counting  distributions, 
namely  to  ir  that  are  (1)  Poisson,  (2)  Binomial  or  (3)  Negative 
Binomial  (Pascal)  (See  also  Sundt  &  Jewell  (1981)  for  genera¬ 
lizations)  . 


The  Binomial  Counting  Distribution 


From  (16)  (17),  we  know  that  for  our  problem  we  want  the 
variance  of  n  to  be  smaller  than  the  mean;  this  suggests  an 
improved  approximation  might  result  from  using  a  Binomial 
counting  density; 


(24) 


T 

n 


rM) 


nn(l-n) 


M-n 


(n=0,l, . . .M) 


with  moments 

(25)  E(h)  =  -M  ;  l/(n)  =  ^M(I-tt) 

For  this  counting  law,  Panjer  (1981)  shows  that  (23)  is 
replaced  by; 


-8- 


g(0)  =  (l-i) M 


(26) 


g(y) 


T  min(y,R) 


x=l 


l (M+l) (x/y) -1 ] f (X) g(y-x) 


(y«l,2, . . .MR) 


so  that  the  recursive  computation  is  still  more  efficient  than 
using  (9)  . 

Note  especially  that  we  are  not  proposing  to  set  M=N,  so  that 
both  parameters  (M,i)  are  available  to  match  (16)  (17).  For  an 
exact  match  of  the  first  two  moments,  we  require  that: 

(27)  M  *  (Iqtmi)2/(Zq^)  j  iM  -  ^ 

(The  reader  may  easily  show  that  MSN) . 

However,  the  Binomial  recursive  algorithm  only  works  for  M 
integer,  so  the  value  obtained  above  must  be  rounded  up  or  down, 
and  then  i  readjusted  to  provide  an  exact  fit  to  the  first  moment. 
The  variance  of  h  and  of  y  will  be  slightly  too  large  (too  small) , 
compared  with  (16)  (17),  if  M  is  adjusted  upwards  (downwards) 

from  the  exact  value.  But  this  error  is  in  general  quite  small 
for  moderate  values  of  N. 

The  Binomial  counting  distribution  also  has  a  good  theoretical 
justification,  for  if  the  original  portfolio  is,  in  fact,  homo¬ 
geneous,  with  and  for  all  i=l,2,...N,  then  we  have 

in  exactly  Binomial,  with  M*N  and  i=qg. 

For  small  inhomogeneities,  if  we  set: 

(28)  m^  »  m  +  ;  q^  =  qg  +  ;  ( i*l , 2 , . . .N) 

where  m  is  defined  in  (14),  and 


q0  =  :q3/N 


(29) 
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wa  find  that,  to  first-order  terms  in  and 

(.30)  t  =  qQ[l  +  (2/N)Z(ui/m)  i  ; 

(31)  M  z  N[1  -  (2/N) Z(ut/m) J  ; 

that  is,  only  the  small  inhomogeneities  in  affect  the 
values  of  M  and  ir. 

If  our  original  portfolio  becomes  quite  large  (N-<”)  , 
but  the  policy  characteristics  (q^m^)  remain  comparable, 
then  (27)  implies  that  M  is  of  order  N  and  will  thus  increase 
without  limit,  but  that  tt  will  remain  relatively  stable. 

This  means  that  we  do  not  expect  that,  in  the  limit,  our 
Binomial  counting  law  will  become  approximately  Poisson 
(which  would  require  M-°°,  t-0,  with  Mit=a)  .  The  justification 
of  the  Poisson  law  thus  requires  other  limiting  conditions. 

Associated  Functions 

As  pointed  out  by  Gerber  (1980) ,  once  a  recursive  procedure 
for  the  density  g(y)  *  Pr(y»y}  has  been  set  up,  it  is  a  trivial 
matter  to  initialize  and  calculate  other  associated  functions. 
The  functions  which  seem  of  most  interest  are  the  complementary 
distribution  function: 

»  y 

Gc(y)  =  Pr(y>y)  =  Z  g(x)  ■  1.0  -  Z  g(x)  , 

x=y+l  x=0 

and  the  "stop-loss  premium11: 

jo  v-1 

.1  ,  (y)  =  Et  (y-y) +  )  -  Z  Gc(x)  =  E(y)  -  I  Gc(x)  . 

si-  x=y  x=0 
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Extension  to  Negative  Discrete  Values 


Aaelson's  and  Panjer's  algorithms  were  developed  only  for 
positive  wi,  which  is  why  the  above  discussion  was  limited  to 
the  sum  of  non-negative  x. .  However,  Sundt  A  Jewell  (1981) 
indicate  how  arbitrary  values,  say  x^  in  the  range  [ -L, . . . , 0 , . . . R ] 
(L,R>0) ,  can,  in  principle,  be  handled  for  Poisson  or  Binomial; 
we  develop  only  the  Binomial  case. 

First  of  all,  (2)  is  replaced  by: 

(32)  fi(x)  *  Pr {x^»x ,  x^O }  , 

and  (26)  is  replaced  by: 

«  T  '.B 

(33)  g (y )  *  -r— :Z  i (M+l) (x/y)  -  l]f(x)g(y-x)  ,  ( -MLSySMR;  y*0) 

~ '  X=*A 

where  A=max ( y-MR, -L) ,  and  B=min(y+ML,R)  . 

g(0)  is  no  longer  calculable  explicitely  from  this  form,  but  both 
g(-ML)  and  g (MR)  are  available  from  first  principles,  and  (33) 
can  be  re-arranged  to  start  the  recursion  at  either  end. 

Starting  from  the  lower  end,  we  obtain: 

( 

1  0  ( y<  -MI.)  and  (y>MR) 

f.  Tf(-L)]M  (y— ML) 

=  1 

I 

LTy+ML) f (-L) j  ~(T  1_1) (y-L)g(y-L) 

C 

+Z  [ (M+l) x-ML-y ] f (x-L) g (y-x)  (otherwise), 
x*l 


(34) 

g(y) 


with  C-min(y+ML,L+R) . 
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Of  course,  if  L  is  very  large,  then  there  are  obvious 
proolems  with  the  accumulation  of  round-off  error,  especially 
if  fl-L)  and  nearby  values  are  small.  One  can  also  imagine 
multi-pass  recursive  procedures,  or  iterative  techniques 
using  (33)  to  resolve  these  numerical-analytic  difficulties. 


Further  Improvements 


One  can  readily  imagine  a  variety  of  further  improvements 
to  the  Binomial  compound  law  to  provide  a  better  approximation: 
for  example,  since  (9)  is  linear  in  the  (it  )  ,  one  could  take 
a  linear  mixture  of  several  counting  distributions,  and  then 
mix  the  results  of  the  corresponding  recursively  calculated 
aggregate  claim  densities;  this  would  enable  matching  higher 
moments  or  other  attributes  of  the  true  density  (4) . 

One  direction  which  we  have  examined  is  to  provide  a  better 
fit  to  the  true  value  of  g(0)=IIp.  ,  which,  as  previously  mentioned, 
is  too  large  in  the  Poisson  case;  the  Binomial  law,  g(0)*^0=(l-,r)  , 
seems  to  give  a  better  numerical  fit,  but  we  cannot  guarantee 
this. 

In  Sundt  &  Jewell  (1981),  it  is  shown  how  to  modify  the 
Panjer  algorithm  so  that  the  new  counting  density  (tM)  can  take 
on  values: 

f  c  +  (l-c)Tfn  (n=0) 

(35)  »i  -  j 

l  (l-o)*n  (n=l, 2 , . . . ) 

V 

where  the  tn  are  Binomial  (w,M).  Alternatively,  one  can  continue 
to  use  (26)  ,  and  mix  the  resulting  density  in  the  obvious  way 
with  the  degenerate  density  at  zero.  This  modified  Binomial 


compound  law  gives  us  three  degrees  of  freedom  (j,t,m). 


Assuming  that  claim  amounts  are  positive,  we  can  match 
g(0)  by 

(36)  g ( 0)  =  Tg  =  •:  +  (l-o)  =  ^Pi 

and  ( 2 S )  becomes 

(37)  £  (h)  =  (l-^)TM  ;  1/  (n)  =*  (1-c)  [  tM(I-t)  +  ] 

These  must  be  matched  numerically  to  the  true  values  (22)  (16)  (17) 

by  iterative  numerical  methods,  which  we  shall  not  describe.  As 
before,  the  integrality  of  M  means  that  we  cannot  exactly  match 
both  the  second  moment  and  g(0) ,  so  that  one  has  to  decide  which 
improvement  is  more  important. 

We  shall  see  in  the  example  to  follow  that  this  modified 
Binomial  provides  only  a  modest  improvement  over  the  Binomial, 
and  suggests  that  further  ref  inements  will  be  of  marginal  value. 


A  Numerical  Example 


To  illustrate  the  effect  of  the  approximation  improvement, 
we  use  a  numerical  example  due  to  Gerber  (1979) ,  in  which  there 
are  N»31  policies,  and  the  random  values  x^  are  either  0  or  a 
"face  value",  c  ^ ,  with  probability  p.^  or  q^,  respectively,  as 
shown  in  Table  X  (The  duplication  of  identical  policies  is 
typical).  Thus  m^-c^  and  v^=0  (i=l , 2 , . . . 31) . 


*1 

Face  Values  c^ 

1 

2 

3 

4 

5 

.03 

2 

3 

1 

2 

- 

.04 

- 

1 

2 

2 

1 

.os 

- 

2 

4 

2 

2 

.06 

- 

2 

2 

2 

1 

Table  I.  Number  of  Policies  with  Indicated  and  c ^ . 
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The  exact  values  of  the  density  g(y) ,  the  complementary 
distribution  Gc(y),  and  the  stop-loss  premium  n  ^(y),  were 
obtained  by  convolving  31  two-point  (0,c^)  densities,  and 
are  given  in  column  three  of  Tables  IV,  V,  and  VI.  From  (6) (7), 
we  find  that  the  first  two  moments  of  the  original  portfolio 
are : 

£(y)  *  4.49  ;  V  (y)  =  15.3003 

and  that  g ( 0 )  =  0.23819 

The  unnormalized  prototypical  claim  density  (13)  used  in 
both  collective  risk  approximations  is  shown  in  Table  II. 


B1SII 

IEI 

2 

3 

4 

5 

1.4  f(x) 

.06 

.35 

.43 

.36 

.20 

Table  II.  Density  of  equivalent  homogeneous  claims. 


The  first  two  moments  of  this  " severity"  density  are: 

ra  *  3.207143  ;  v  *  1.207092  . 

Thus,  from  (16)  (17),  the  "counting"  density  moments  for 
an  exact  fit  of  a  collective  risk  approximation  must  be: 

£ (n)  =  1.4  ;  V(n)  =  1.323224 

Three  approximations  were  computed  using  recursions  (23)  and 
(26)  and  the  method  of  (35), giving  the  numerical  matching  shown 
in  Table  III. 


Exact 

Values 

— 

Approximations 

Poisson 

Binomial 

Modified 

Binomial 

E(y) 

4.49 

4.49 

4.49 

4.49 

V  (y) 

15.3003 

16.0900 

15.3146 

15.3COJ 

Pr(y»0}*g(0) 

0.23819 

0.24660 

0.23714 

0.23809 

Table  III.  Value  Matching  for  Numerical  Example. 
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In  the  Poisson  approximation,  -=1.4  fixed  E(y)=4.49  as 
desired  but  1/  (y)  *16 . 0900  and  g(0)=0. 24660  are  significantly 
too  large.  Results  using  the  recursion  (23)  (part  of  which 
were  also  given  in  Gerber  (1979))  are  shown  in  column  two 
of  Tables  IV,  V,  and  VI. 

For  the  Binomial  counting  distribution,  an  exact  match  of 

the  first  two  moments  would  require  M=25. 528480  and  t=0. 0548400. 

Rounding  up,  we  select  integer  M=26,  and  adjust  *=0. 0538462  to 

keep  E(y)=4.49.  1/ (y)  =15 . 3146  is  still  significantly  close  to 

the  exact  value  of  15.3003,  but  g(0)=0. 23714  is  now  less  than 

the  true  value.  Note  that  the  range  of  the  Binomial  approximation 

extends,  in  principle,  to  5x26=130,  whereas  the  largest  possible 

total  claim  sum  of  the  original  portfolio  is  only  97.  However, 

reference  to  Table  V  shows  that  the  probability  of  a  claim 

-9 

larger  than  40  is  already  of  order  10  I 

For  the  modified  Binomial  approximation,  we  must  use  (36) (37) 
to  find  the  parameter  values  to  match  the  first  two  moments  and 
g(0) ;  These  turn  out  to  be  M=21. 737130,  t=0 . 0648672 ,  and 
o=0 . 00711084 .  Rounding  up,  we  set  M=22,  and  readjust  the  other 
values  to  match  the  mean  and  variance,  giving  finally  -=0.064055 
and  j=0. 00653874 .  As  can  be  seen  from  Table  III,  the  resulting 
mismatch  in  g(0)  is  quite  small. 
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g(y)  =  Priy=y; 


APPROXIMATIONS 


B 

RESULT 

POISSON 

BINOMIAL 

i 

MODIFIED 

BINOMIAL 

0 

0.23819 

0.24660 

0.23714 

0.23809 

1 

0.01473 

0.01480. 

0.01504 

0.01494 

2 

0.03773 

C.  08675 

0.08818 

0.08762 

3 

0. 11318 

0.11122 _ 

0.11313 _ 

0.11246 

4 

0.11071 

0.11040 _ 

0.11256 _ 

0.11206 

5 

0.09633 

0.09286 

0.09507 

0.09492 

6 

0.06155 

0.06101 _ 

0.06315 

7 

0.06902 

0.06543 

0.06732 

0.06759 

3 

0.05482 

0.05458 

0.05589 _ 

0.05613 

9 

0.04315 

0.04132 _ 

0.04197 

0.04217 

10 

0.03011 

0.03058 

0 . 0307_1 

0.03086 

q 

0.02353 

0.02331 

0.02311 _ 

0.02321 

0.01828 

0.01834 

0.01797 

0.01802 

m 

0.01251 

0.01315 

0.01265 

0.01266 

14 

0.00871 

0.00922 

0.00866 

0.00865 

15 

0.00591 

0.00650 

0.00596 

0.00593 

16 

0.00415 

0.00460 _ 

0.00411 

0.00408 

17 

0.00272 

0.00318 _ 

0.00277 

0.00273 

13 

0.00174 

0.00212 _ 

0.00179 

0.00176 

19 

0.00112 

0.00141 _ 

0.00115 

0.00112 

20 

0.00071 

0.0009£ 

0.00073 

0.00071 

30 

8.63294X10-6 

3 . 98500xl0~6 

3.51483xl0-6 

40 

Table  IV.  Total  Sum  Densities  in  Example  (Differing  digits 
underlined) . 
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c 

G  (y)  *  Pr  ;y>y 


APPROXIMATIONS 


y 

RESULT 

POISSON 

BINOMIAL 

MODIFIED 

BINOMIAL 

0 

0.76181 

0.75340 

0.76286 

0 . 7619_1 

1 

0.74707 

0.73861 

0.74782 

0.74696 

2 

0.65934 

0.65135 

0 . 6  596_4 

0.65934 

3 

0.54615 

0.5406J _ 

0.54651 

0.54688 

4 

0.43544 

0.43395 

0.43482 

5 

0.33912 

0.33888 

0.33990 

6 

0.27757 

0.27637 

0.27597 

0.27675 _ 

7 

0.20855 

0.21094 

0.20865 _ 

0.20916 

8 

0.15373 

0 . 156_36 _ 

0.15276 

0.15303 _ 

9 

0.11058 

0.11504 

0.11079 

0.11086 

0.08048 

0.08446 

0.08008 

0.08000 

H 

0.05695 

0.06115 _ 

0.05696 

0.05679 

H 

0.03866 

0.04281 

0.03899 

0.03877 

B 

0.02615 

0.02966 

0.02635 

0.02611 

14 

0.01744 

0.02044 

0.01769 

0.01746 

15 

0.01153 

0.01394 

0.01173 

0-01153 

16 

0.00738 

0.00934 

0 . 0076_2 

0.00745 

17 

0.00467 

0.00617 

0.00485 

0.00472 _ 

18 

0.00292 

0.00404 

0.00306 _ 

0.00296 

19 

0.00181 

0.00263 

0.00192 

0.00184 

20 

0.00110 

0.00169 

0.00118 

0.00112 

" 

1  30  1 

3.49840xl0-6 

12 . 46  21xl0~6 

4 . 87524xl0~6 

4.16710X10'6 

J  40 

3.10833xl0-9 

45.5298X10*9 

7.42S41X10'9 

5 . 2601 3xl0~9 

Table  V.  Complementary  Distributions  in  Example  (Differing  digits 
underlined) . 
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:sl(y) 

*  El  (y-y)i’i 

EXACT 

APPROXIMATIONS 

D 

RESULT 

POISSON 

BINOMIAL 

MODIFIED 

BINOMIAL 

4.49000 

4.49000 

4.49000 

4.49000 

3.72819 

3.73660 

3.72714 

3.72809 

2.98112 

2.99799 

2.97932 

2.98113 

2.32179 

2.34614 

2.31968 

2.32179 

1.77563 

1.80551 

1.77317 

1.77491 

5 

1.34019 

1.37527 

1.33922 

1.34009 

6 

1.00106 

1.03790 

1.00034 

1.00019 

7 

0.72350 

0.76153 

0.72437 

0.72345 

8 

0.51495 

0.55059 

0.51572 

0.51428 

9 

0.36122 

0.39423 

0.36296 

0.36125 

0.25064 

0.27919 

0.25217 

0.25039 

Hi 

0.17017 

0.19472 

0.17209 

0.17039 

D 

0.11322 

0.13357 

0.11513 

0.11360 

H 

0.07456 

0.09076 

0.07614 

0.07483 

14 

0.04840 

0.06110 

0.04979 

0.04872, 

15 

0. 03096 

0.04065 

0.03210 

0.03126 

16 

0.01943 

0.02671 

0.02037 

0.01973 

17 

0.01205 

0.01737 

0.01276 

0.01228 

18 

0.00738 

0.01120 

0.00791 

0.00756, 

19 

0.00446 

0.00716 

0.00485 

0.00460 

20 

0.00265 

0.00453 

0.00293 

0.00276 

30 

|  29.7953xl0-6 

8.88376xl0-6 

40 

5.72441xl0"9 

101. 020xl0~9 

14 . 6686xl0-9 

Table  VI.  Stop-loss  Premiums  in  Example  (Differing  digits 
underlined) . 
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Comparison  of  the  different  results  for  the  density,  g(y), 
in  Table  IV  shows,  as  expected,  that  none  of  the  approximations 
is  a  particularly  good  point  estimator;  because  of  the  differences 
in  the  models,  the  approximations  are  forced  to  fluctuate  above 
and  below  the  exact  density.  The  modified  Binomial  is  generally 
better  than  the  Binomial,  which  is  generally  better  than  the 
Poisson,  although  this  is  by  no  means  uniformly  true. 

However,  when  we  examine  the  complementary  distributions, 

Gc(y),  in  Table  V,  the  approximations  become  more  stable,  and 
the  Binomial  is  always  better  than  the  Poisson,  except  for  y=6 . 

The  modified  Binomial  is  uniformly  best  only  from  y*12  onwards. 

The  approximations  to  the  stop-loss  premiums,  :isl(y),  in 
Table  VI,  are  even  more  stable,  and  show  clearly  the  value  of 
matching  the  second  moment  for  this  "tail  of  the  tail".  The 
Poisson  is  always  worst,  and  the  modified  Binomial  always  best, 
except  at  y=6. 

These  remarks  can  be  more  easily  visualized  in  Figures  1,  2, 
and  3 ,  which  show  the  percentage  error  in  each  approximation  for 
the  functions  of  interest.  In  addition  to  the  remarks  above, 
it  is  of  interest  to  observe  the  inevitable  degradation  of  all 
approximations  at  large  values  of  y.  It  can  be  shown  theoretically 
(Biihlmann,  et  al.  1977)  that  the  Poisson  approximation  gives 
too  conservative  (large)  a  value  for  the  stop-loss  premium  for 
all  values  of  y.  Our  example  suggests  all  of  these  approximations 
are  eventually  "too  dangerous"  in  the  tails.  However,  it  should 
be  remembered  that  the  actual  values  of  the  probabilities  and 
of  the  absolute  errors  are  quite  small  above  y*20. 
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(54.2%) 


FIGURE  2 

PERCENTAGE  ERROR  IN  APPROXIMATIONS 
TO  COMPLEMENTARY  DISTRIBUTION  G  (y)  VERSUS  y 


-2 


Conclusions 


Naturally,  only  limited  conclusions  can  be  drawn  from  a 
single  computational  example.  However,  we  believe  that  the 
Binomial  compound  law  is  a  significantly  better  approximation 
to  the  distribution  of  the  original  heterogeneous  portfolio 
than  the  traditional  Poisson  compound  law  approximation; 
furthermore,  it  can  also  be  computed  recursively  with  little 
increase  in  difficulty.  There  also  seems  to  be  evidence  that 
the  slight  additional  work  to  set  up  the  modified  Binomial 
compound  law  approximation  will  be  worthwhile  if  more  accurate 
values  of  the  complementary  distribution  or  the  stop-loss  premium 
are  desired  in  the  tails. 

It  remains  to  be  seen  whether  there  are  significant  differences 
between  these  approximations  for  real  risk  portfolios,  where  N 
and  R  are  both  large,  and  where  round-off  error  accumulation  may 
become  important  in  any  recursive  method.  There  have  been  some 
claims  that  other  approximation  methods  or  fast  Fourier  transforms 
may  De  competitive  under  these  conditions. 

Finally,  we  must  keep  in  mind  the  ever-increasing  capabilities 
of  digital  computers,  and  the  fact  that  many  real  portfolio 
distributions  can  best  be  calculated  directly. 
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